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1. INTRODUCTION 

There has been a lot of interest in neutrino interac- 
tions in the few GeV energy region, coming from neu- 
trino oscillation experiments and their need to better 
constrain systematic errors. Despite a vivid activity in 
the field [2, Q the neutrino-nucleon cross sections un- 
certainties are still around 20-30%. There are several 
reasons for that: 

(i) in all the cross section measurement experiments the 
neutrino flux was known with poor precision, both the 
overall normalization and the spectrum; 

(ii) in all recent attempts to measure neutrino cross sec- 
tions, interactions occur on nuclei with many compli- 
cations coming from nuclear effects. 

For neutrino energies around 1 GeV the most 
abundant reaction is the charged current quasi-elastic 
(CCQE) scattering: vi + n ^ + p, where I e {e,^}. 
It is the most important process in the investigation of 
the oscillation phenomenon, e.g. in the T2K experiment 
After applying G-symmetry arguments the weak 
transition matrix element of the hadronic current enter- 
ing the CCQE cross section formula can be expressed 
in terms of two vector and two axial form-factors. The 
former ones are known from electron scattering exper- 
iments due to the CVC (conserved vector current) hy- 
pothesis. The latter are proportional to each other due 
to the PCAC (partially conserved axial current) hypoth- 
esis, hence only one independent axial form factor re- 
mains. It is usually assumed to have a dipole form (Eq. 
IS]) parametrized by the axial mass Ma- 

Recent neutrino cross section measurements, notably 
the high statistics muon double differential cross sec- 
tion results from the MiniBooNE (MB) experiment [1], 
suggest significantly larger Ma values than estimations 



based on older deuterium target neutrino measurements 
or on the pion electro-production data Q 

In the MB experiment the signal events (CCQE-like) 
were those with no pion in the final state. CCQE sample 
of events was obtained after subtraction of an estimated 
(based on the data / Monte Carlo (MC) comparison) 
contribution from pion production and absorption. 

In the last years it is becoming clear that the MB 
analysis neglected a large two body current contribu- 
tion to the cross section [3, @] giving rise to events that 
can be easily confused with genuine CCQE events un- 
less one carefully investigates the final state nucleons. 
Older and recent Ma measurements can be consistent, 
because in the case of neutrino-deuteron scattering this 
contribution is small [9]. 

Several groups attempted to explain the MB CCQE 
double differential cross section data with models con- 
taining a significant contribution from the np—nh mech- 
anism (n particles and n holes, n > 2, i.e. n nucleons in 
the final state). The np — nh mechanism is also called 
meson exchange current (MEC), multi-nucleon knock- 
out or two-body current contribution. The IFIC group 
performed a fit to the MB CCQE data and obtained 
Ma = 1.077± 0.027 GeV ^10,]. A good qualitative agree- 
ment was reached within the Lyon group model 11[ as 
well as the optical potential model 11 21 and slightly worse 
within the superscaling approach Il3l and also the trans- 
verse enhancement (TE) model [ij, [l3| ■ Even though 
the theoretical models are able to explain the large Ma 
result of the MB measurement, but their predictions for 
the size of the np — nh contribution in some cases (2D 
bins as used by the MB collaboration) differ by a factor 
of two. 

The theoretical models of the MEC contribution give 
quite different estimates of a significance of the effect in 
the case of antineutrino scattering. Recently the MB col- 



2 



laboration published the first antineutrino large statis- 
tics CCQE cross section results UM. The data has been 



already analyzed by IFIC group [17[ . In [18| the ratio of 
CCQE-like cross sections (defined as explained above) 
for neutrinos and antineutrinos was discussed as a func- 
tion of neutrino energy. This data can allow to discrimi- 
nate between the models^ . The results seem to favor the 
IFIC model but the data suffers from large errors and 
only the superscaling approach seems to be in trouble. 

The importance of the CCQE reaction gives a strong 
motivation to look for alternative ways to evaluate Ma 
and/or investigate the size of the two body current con- 
tribution. An interesting option is provided by the neu- 
tral current elastic (NCEL) reaction: vi + N ^ vi + N , 
where N denotes proton or neutron. 

As explained in Sect. [2] the basic theoretical frame- 
work to investigate NCEL scattering is similar to the 
one used in CCQE. The hadronic current is expressed 
in terms of vector and axial form factors. An inter- 
esting feature is that the NCEL form factors are linear 
combinations of the form factors present in the CCQE 
reaction with the addition of terms sensitive to a strange 
quark content of nucleons, see Eq. [B] and [T] Thus, the 
NCEL scattering data allows to estimate both Ma and 
the strange quark contribution to the form factors. In 
fact, most of the interest in NCEL reaction comes from 
its potential to measure strangeness of the nucleon. 

Historically, the first estimations of the nucleon 
strange quark content were based on the analysis of the 
deep inelastic scattering (DIS) data from CERN, SLAC, 
DESY, and Jefferson Lab, for a review see [ISl- There is 
an interesting relation between the quantities measured 
in DIS and the nucleon form factors. For instance, from 
the asymmetry measurement of longitudinally polarized 
leptons off proton one can establish the first moment 
of the proton structure function gi. Assuming a naive 
parton model and taking into consideration the isospin 
symmetry as well as an appropriate sum rule, it can be 
shown that the fraction As of strange quarks and anti- 
quarks which contribute to the total spin of the proton, 
corresponds to gf = G\{Q), see Eq. [TTJ Old mea- 



surements suggested the nonzero value of As, see |2l| . 
however, more recent estimations give results consistent 
with zero 

The axial strange form factor can be determined from 
the NCEL neutrino-nucleon or neutrino-nucleus scatter- 
ing data. The first results were obtained by the BNL 
E734 experiment [i^]- The analysis of the vp — vp 
and vp -> vp data gave: Ma = 1-06 ± 0.05 GeV and 
= —0.15 -I- ±0.09. For later discussions of these data 



see 



24, 2 



^ Notice, however, a criticism of the adopted energy unfolding 
procedure explained in [Tgll . 



About three years ago the MB collaboration measured 
the flux averaged NCEL differential cross section in 
on the CH2 target The MB proposed also an ob- 
servable, a proton enriched sample of events, sensitive 
to the strange axial form factor. The first measurement 
relied on the MB detector ability to analyze the scin- 
tillation light even in the absence of Cherenkov light 
from the final state muon. The dependent differ- 
ential cross section on an almost isoscalar target de- 
pends very weakly on 17^, so a conclusive simultaneous 
extraction of both Ma and gf parameters is impossible. 
In the second measurement, a high energy (above the 
Cherenkov radiation threshold) proton enriched sample 
of events were selected. The MB performed two sep- 
arate parameter extractions. Assuming 5^ = the 
value of M"/^ = 1.39 ± 0.11 GeV was obtained. Taking 
M^J^ = 1.35 GeV from the CCQE analysis [l lead to 
the estimation of g\ = 0.08 ± 0.26. 

In the MB analysis the values of Ma used in model- 
ing scattering off carbon (the value is floating) and off 
free protons, were different. In the case of proton the 
Ma was fixed at 1.13 GeV. For carbon the axial mass 
was treated as an effective parameter. A large value of 
effective Ma is expected to account for the MEC con- 
tribution present in neutrino-carbon scattering but ab- 
sent in neutrino-proton scattering. MB considered also 
a non-zero value of the k parameter enhancing the Pauli 
blocking effect at low Q^. Because the CCQE data are 
consistent with k = and the physical interpretation of 
this parameter is unclear we did not include it in our 
analysis. 

The data from [1] was already discussed in several 
phenomenological papers. Butkevich and Perevalov 26 1 
investigated an impact of using a more realistic nu- 
clear model in the MB data analysis. The applica- 
tion of a relativistic distorted wave impulse approxi- 
mation (RDWIA) jl] model lead to the values Ma = 
1.28 ± 0.05 GeV and g^ = -0.11 ± 0.36, consistent with 
those reported in [ij. The spectral function formalism 
was applied by Ankowski in ^2S] . with a result that the 
shape of the MB measured distribution of events in 
is reproduced with M'^^ = 1.23 GeV. However, there 
is a 20% discrepancy for the overall normalization with 
the MB results (the measured cross sections are larger) . 
Meucci, Gusti and Pacati analyzed the predictions of 
four nuclear models It turned out that relativistic 
Green function (RGF) model is able to reproduce MB 
NCEL data with a value of Ma close to those obtained 
in deuterium-target experiments. Apparently, the opti- 
cal potential discussed in the RGF formalism accounts 
for the MEC contribution. 

The main purpose of this paper is to analyze the MB 
NCEL data with the NuWro Monte Carlo event gener- 
ator [sol with the MEC contribution and final state in- 
teractions (FSI) effects implemented there. Out of three 
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MEC models available in NuWro, in the current analysis 
we decided to apply the (TE) model [14|. We adopted 
the same approach as proposed in the original paper for 
the charged current scattering: the MEC contribution is 
defined by the same modification of the magnetic form 
factor [31[. The existing models for the MEC contribu- 
tion give predictions for the final muons only. To obtain 
final state nucleons, we use a procedure proposed in jsil ]. 
This approach is universal and was applied also to two 
other MEC models implemented in NuWro. We expect 
that the differences between the three MEC models pre- 
dictions for nucleons are small and we decided to work 
with the TE model because of its simplicity. 

Another important difference with respect to previous 
studies of the MB NCEL data lies in a treatment of 
nuclear effects. In our analysis we follow closely the 
approach proposed by Perevalov in his PhD thesis 33 1. 
We distinguish three scenarios for nucleons arising from 
the initial NCEL interaction on carbon: 

(a) a proton leaves nucleus without reinteractions; 

(b) a proton from the primary interaction is subject 
to reinteractions and not specified nucleon/nucleons 
are knocked out; 

(c) NCEL interaction occur on a neutron. 

We also take into account NCEL scattering on hydro- 
gen and the irreducible background (a pion production 
in primary vertex, which is absorbed during ESI). In all 
five cases we calculate visible energy using response ma- 
trices provided in |33| . The probabilities of scenarios (a) 
and (b) depend on the details of the cascade models im- 
plemented in NUANCE (used by the MB Collaboration) 
and in NuWro. In the MEC events there are typically 
two nucleons after a primary interaction and both prop- 
agate independently through nucleus. For those events 
we had to invent a procedure connecting their kinetic 
energies with the visible energy (see Sect. 01). 

Our main result is the conclusion that after taking 
into account the MEC contribution, the MB NCEL ex- 
perimental data can be reproduced with the Ma value 
that is very close to those obtained in old deuterium 
target experiments. The procedures estimating the vis- 
ible energy was verified by repeating the MB analysis 
without the MEC contribution. Despite the difference 
in modeling of the FSI effects by NUANCE and NuWro, 
our Ma estimation is similar to the value obtained by 
the MB collaboration. 

Our paper is organized as follows: in Sec. [Da general 
description of our theoretical model for the NCEL reac- 



tion is given; in Sec. |3]the main features of the NuWro 
generator are summarized; Sec. |4] contains a details of 
the data analysis and an energy unfolding procedure; in 
Sec [5] we present main results and our conclusions can 
be found in Sec. [51 

2. ELASTIC NEUTRAL CURRENT REACTION 
FORMALISM 

Let us consider neutral or charged current neutrino- 
nucleon scattering: 



iyik)+Nip)^l\k') + N'{p'), 



(1) 



where N, N', I' denote the initial and final nucleon and 
the final lepton with four momenta: p, p' and k' , re- 
spectively. The transfer of four momentum is given by 
q^' = k^'- k'" = (cj, q), = ^q^. 

In the Born approximation the scattering matrix ele- 
ment factorizes into leptonic and hadronic parts: 



*^cc,nc 



Q'^M^,Ml (2) 



where fee = cos 9c, fnc = 1, Oc is Cabibbo angle, G_f is 
the Fermi constant, while and /ijf^ „^ are the expecta- 
tion values of the leptonic and hadronic currents. The 
leptonic part is easy to derive: 



= u(fc')7^(l - 75)u(fc), 



(3) 



while getting the hadronic contribution requires an extra 
phenomenological input. The most general form of the 
expectation value of the hadronic current reads 



h^{q) = u{p')r>^{q)u{p), 



(4) 



where F'' is an effective hadronic vertex modeled by form 
factors. 

In order to construct the neutral current hadronic 
vertex for the elastic neutrino-nucleon scattering, one 
has to follow a general pattern given by the Standard 
Model. Indeed, if one takes into consideration u, d, and 
s quarks then one can construct the one body quark 
current, which accounts also for the strangeness contri- 
bution. The CVC theory, the PCAC hypothesis and the 
SU (2) isospin symmetry relate the neutral current form 
factors to those present in electormagnetic and charged 
current hadronic matrix elements 20[ . The neutral cur- 
rent hadronic vertex reads 



J 



NC,p(n) 



NC,p{n) 



2M 



(5) 
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where indices p and n refer to proton and neutron. The neutral current form factors can be expressed as: 

F,7.^(")(g2) = ±\{Fl,{Q')-Fl,{Q')]-2,in^ewFl^-\Q')-\Fl,{Q% (6) 
g^c.p(n)^Q2) ^ ±Iga{Q') - \g\{Q% (7) 



+/— signs refer to proton/neutron, 9w is Weinberg an- 
gle, sin^ 9w = 0.23120. 2""* ^""^ proton (neutron) elec- 
tromagnetic form factors, Ga is the axial nucleon form 
factor, 

9A 



Ga{Q') = 



1 



Ml 



9A = 1.267. 



(8) 



Ff 2 and are the vector and the axial strange form 
factors. 

The electromagnetic form factors are obtained from 
the analysis of the elastic eN scattering data (for a re- 
view see (15). They can be expressed by the electric 
and magnetic form factors: 

Ff(")(g2) = 



4M^ 



Q2 + 



AM' 



(9) 

Gt\Q')^Gt\Q')\m 



where M = \ {Mp + M„) is an average nucleon mass. 
In the Breit frame, at low Q^, the electric (G^^ 



and magnetic (G^^"-*) nucleon form factors are related 
to the electric and magnetic charge distributions inside 
the nucleon. For instance, for the electric proton form 



factor G^(g2) 

,dGl 



6 



dQ2 



G^(0) + %^Q^ + O(g4), where (r^) = 
is the mean-square radius of the charge 



distribution. 

Analogically, we can introduce the electric, G|;, and 
magnetic, G^, isoscalar strange form factors. In the 
first approximation one can assume that "the effective 
size" of the proton and neutron strange sea is the same, 
therefore we use the same strange form factors. The 
latest global analysis of the elastic PV ep and vN BNL 
scattering data suggests that the vector strange form 
factors are consistent with zero [2^ . consequently, in 



our analysis we set G|; 



0. 



In order to estimate the axial strange form factor we 
assume that the radius of the strange sea is comparable 
with the axial "charge" radius of the proton. Therefore 
we take the same dipole form of the strange axial form 
factor as in Eq. [51 



Gi(g2 



9a 



(11) 



and only g\ has to be established from the data 

In our numerical analysis we use BBBA05 vector form 
factors '35']. However, for completeness, we also checked 
a possible impact of the form factors corrected by two- 
photon exchange effect [3£ 
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NUWRO MONTE CARLO EVENT 
GENERATOR 



3.1. Generalities 

NuWro is a Monte Carlo event generator developed at 
the Wroclaw University. It simulates neutrino-nucleon 
and neutrino-nucleus charged and neutral current in- 
teractions including: (quasi-) elastic scattering, pion 
production through A resonance (with some contribu- 
tion from a non- resonant background), more inelas- 
tic processes (in neutrino MC community often called 
deep inelastic scattering) and coherent pion produc- 
tion. NuWro covers neutrino energies from ~ 100 MeV 
to TeV. There are two nucleus models implemented in 
NuWro: the Fermi gas (FG) model and an approach 
using spectral function. All hadrons from primary ver- 
tex are propagated through the nuclear matter using the 
NuWro cascade model. A detailed description of NuWro 
and its physical models can be found in [3(| . NuWro is 
an open-source project and the code is available from 
the repository [38|. 



3.2. Meson Exchange Current Contribution 

There are three MEC models available in NuWro: the 
IFIC model, the Lyon group model and the TE model. 
In each of them the double differential cross section for 
the final state muon is the main external input and the 
hadronic part is modeled using the scheme described 
below. 

In our analysis we use the transverse enhancement 
model described in [iJl- Initially, the TE model was 
introduced to describe CCQE interactions, but it can 
be applied to neutral current reactions in a similar way 
The MEC contribution to the NCEL scattering 
cross section on carbon is introduced by a modification 
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FIG. 1: [Color online] Total cross section per nucleon as 
a function of a neutrino energy for NCEL scattering (with 
Ma = 1350 MeV) and TE model on cabon. 



of the vector magnetic form factors: 




exp 



(12) 

where A = Q GeV"^ and B = 0.34 GeV^. Using G^f in 
the NCEL differential cross section formula one obtains 
a cross sections for NCEL scattering with the MEC con- 
tribution included. The MEC cross section is calculated 
from the difference between cross sections with modified 
and standard magnetic form factors: 



2„MEC 



dqduj 



1 

2 dqduj 

^2^NCE ^ 



dqdjjj 



M. 



dqduj 



j2 NCE 



.(13) 



Note that the value of the axial mass is set to be 
M^EC ^ 2^014 MeV, as assumed in Q. The total 



cross section predictions for the model are shown in Fig. 
[TJ At typical MB flux neutrino energies, 700 MeV, 

the MEC contribution is about 15% of the NCEL cross 
section. 

For a description of nucleons a multinucleon ejection 
model from 32| is used. Correlations between two initial 
nucleons are neglected, which significantly simplifies the 
calculations and does not have noticeable infiuence on 
the final nucleons kinematics. With this assumption our 
procedure is the following: 

1. set randomly the four-momenta {pi and P2) of the 
initial nucleons from the Fermi sphere with radius de- 
termined by the local nuclear density; 

2. calculate the four- momentum of the hadronic system 

W = pi + P2 + q, 



where q is the four-momentum transfered to the 
hadronic system; 

3. repeat steps 1. and 2. until the invariant mass in the 
hadronic center of mass system is larger than mass of 
two nucleons {W'^ > 4M^); 

4. make Lorentz boost to the hadronic center of mass 
system; 

5. select isotropically momenta of two final nucleons; 

6. boost back to the laboratory frame. 

Each of final nucleons undergoes final state interac- 
tions. 

In the case of neutrino charged current MEC reactions 
there may be neutron-neutron or proton-neutron pair 
in the initial state. The probability of the mixed pair 
is given by a free parameter pcc, with default value 
Pec = 0.6. 

For neutral current MEC interactions each isospin 
pair is possible (n-p, n-n, p-p). To keep the same pro- 
portion between n-n and p-n pairs we introduced the 
parameter patc = (2/pcc ~ 1)~^ giving the probabilty 
of the mixed pair to be present in a neutral current MEC 
reaction. Both n-n and p-p pairs are assumed to be 
equally likely. 

We checked that changes to the parameter pcc have 
negligible impact on the final results. 



3.3. Nucleus Model 

In order to describe the carbon nucleus, present in 
CH2 molecule, the local relativistic Fermi gas model is 
applied. The impulse approximation is assumed for all 
dynamical channels but the MEC and coherent ones. 
Final state nucleons present in the primary vertex are 
propagated through nuclear matter as modeled by the 
NuWro cascade: 

1. the nucleons are assumed to be in the potential well 
of depth V = Vq -\- Ep, where Ep is the Fermi energy 
and Vq = 7 MeV; 



2. a formation zone is applied, see [3C 
the following steps are repeated: 



for details, and 



3. nucleon's free path (A) is drawn from the exponential 
distribution taking into account nucleon- nucleon cross 
section and the local nuclear density; 

4. if A < Xmax = 0.2 fm the nucleon is propagated by 
A, the interaction kinematics is generated and a check 
for Pauli blocking occurs to decide if the interaction 
happened; 

5. if A > Xmax the nucleon is propagated by Xmax- 
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In the cascade the Fermi motion of target nucleons 
is taken into account, so each nucleon from the nucleus 
brings an extra kinetic energy. When the nucleon leaves 
the nucleus, V is subtracted from its kinetic energy, set- 
ting it on-shell. If its kinetic energy is smaller than V, 
the nucleon is stuck inside the nucleus. 

For more detailed description of the NuWro cascade 
model see |30|. 



4. DATA ANALYSIS 

We are going to discuss two samples of data provided 
by the MB. The first one is called the NCEL sample and 
contains the total reconstructed kinetic energy distribu- 
tion of all nucleons in the final state normalized to the 
number of events seen in the detector. 

The second sample of data, the high energy NCEL 
sample, describes numerator and denominator of the ra- 
tio: 

f] = X{iyp^ iyp)/X{iyN ^ lyN), (14) 

where X denotes a contribution from a special class of 
events, called single proton or proton enriched. Those 
are the events with visible Cherenkov light and kinemat- 
ical cut on the proton angle < 60°. The biggest contri- 
bution to those events comes from the NCEL scattering 
on protons, which then do not undergo reinteractions. 
In the case of multiple proton events, the energy of indi- 
vidual proton is in general too low to produce Cherenkov 
light. Even if a high energy proton appears in the mul- 
tiple proton event, it has usually larger scattering angle 
than proton unaffected by FSI. The denominator (X) 
denotes a contribution from all NCEL interactions. 

Both samples are presented as a function of recon- 
structed energy^ (v), that is seen in the detector. To 
compare those data with the theoretical predictions 
given in terms of true kinetic energy (/i), one needs an 
unfolding procedure, allowing a passage from ^ to v. 

In our analysis, we follow the main steps of the proce- 
dure adopted by the MB collaboration as it is described 
in [SS*]. We had to propose a consistent treatment of 
the MFC events, which were not considered in the MB 
analysis. In the next subsections we describe the orig- 
inal MB unfolding procedure and the one used in this 
paper. We checked, that in the case of the absence of 
MEC, both procedures lead to similar results. 



^ We keep the original notation from [3^ . 



4.1. MiniBooNE Procedure 

The energy unfolding procedure is the same for both 
samples of data. Five different types of signal giving 
contribution to the final distribution are distinguished: 

1. NCEL on hydrogen; 

2. NCEL on a proton from carbon unaffected by FSI 
effects; 

3. NCEL on a proton from carbon with FSI effects; 

4. NCEL on a neutron from carbon; 

5. irreducible background (pious produced in a primary 
vertex and absorbed during FSI). 

For each type of the signal events, fc = 1, 2, 5, there 
is a proper response matrix (i?^*^-*), which simulates the 
energy smearing, the detector efficiency and defines the 
relation between true and reconstructed energy distri- 
butions: 

i 

i?^'^) are either 51 x 51 or 30 x 30 matrices for NCEL 
sample and NCEL high energy sample, respectively. The 
columns of matrices label the true kinetic energy bins 
and rows label the reconstructed energy. There are 50 
bins starting from MeV to 900 MeV plus an extra 
overflow bin for the NCEL sample in the true kinetic 
energy. For NCEL high energy sample there are 28 bins, 
starting from SOOMeV to 900MeV plus the underflow 
and the overflow bins. 

To obtain the reconstructed kinetic energy distribu- 
tion and compare it with the data from fj] one should 
follow the following steps: 

1. use a theoretical model and calculate the flux- 
averaged distributions for five different types of signal 
events using the same bins as in the response matrices; 

2. use the proper response matrices to translate each his- 
togram to the reconstructed kinetic energy distribu- 
tion; 

3. sum all the histograms and add the background events 
(^-^BKG contains dirt, beam-unrelated, and other back- 
grounds) provided by the MB to get the total recon- 
structed energy spectrum: 

+ E^gVf^+E^iM'^ (16) 

i i 

+ E^svr^+-f"^ 
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4. use the provided error matrices (M^ ) to calculate x^: 



X 



Ml 



^ATA _ ^C^ 



(17) 

Note that unlike in the CCQE MB data published in 
here the normalization error is included in the error 
matrix M,,-. 



4.2. Our Procedure 

An alternative way to convert the true kinetic energy 
to the reconstructed one is to translate it on the event 
by event basis. Each value of the true kinetic energy 
is related to the corresponding column in the response 
matrix. This column gives a distribution with an in- 
formation how given true energy is smeared out in the 
detector, normalized to the efficiency. To obtain recon- 
structed kinetic energy distribution one proceeds as fol- 
lows: 

1. for each event calculate total true kinetic energies of 
all nucleons in the final state (/i) and get a bin number 

2. find the kind of signal (fc), see Sect. 14.11 

3. choose j-th column of the BP^^ response matrix as the 
probability distribution; 

4. use the MC method to decide if the event is accepted 
(according to efficiency) and, if so, what energy would 
be visible in the detector. 

This strategy allows us to treat each nucleon from 
MEC event individually, which is required in the NuWro 
based analysis of the data. 



MEC Events m NCEL Sample 
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FIG. 2: [Color online] Number of events as a function of the 
reconstructed kinetic energy of all nucleons in a final state 
measured by MiniBooNE together with our best fits for two 
cases: the MEC contribution was not included (solid line) 
and the MEC contribution was included (dashed line) 



neutrons in primary vertex or for the third signal (NCEL 
on proton from carbon affected by FSI effects) in other 
cases p9i] . 



MEC Events in High Energy NCEL Sample 

In a MEC event the energy transfer is shared by two 
nucleons, so the probability that there will be a proton 
with energy high enough to produce Cherenkov light is 
very low. If there is a proton with the azimuth angle 
below 60°, the second response matrix (NCEL proton 
from carbon unaffected by FSI effects) is applied to its 
true kinetic energy to check if it is over the Cherenkov 
threshold and decide if it contributes to the numerator 
of the ratio 77. 

In both, the numerator and the denominator, we use 
the reconstructed energy of the event similarly as for the 
NCEL sample. 



In the MB analysis there are no MEC events included 
and no response matrices were prepared for them. To 
take MEC events into account, we had to express them 
in terms of five signals defined by MB. 

Naive interpretation may suggest a treatment of each 
nucleon from the MEC events separately. However, it 
would be incorrect, because any of two individual nu- 
cleons may not create enough PMT hits to be visible in 
the detector, but together they can make it (sof . 

In our analysis we treat both nucleons from a MEC 
event together and sum up the energies of all nucleons in 
the final state as if they were coming from one nucleon. 
After that, we apply response matrix for the fourth sig- 
nal (NCEL on neutron from carbon) if there were two 



5. RESULTS 

The defined in Eq. [T7]was minimized for two mod- 
els: with and without the MEC contribution. The pro- 
cedure described in Sec. 14.21 allowed us to treat the nu- 
cleons from MEC events properly, however, we checked 
that the results (for the case without the MEC contribu- 
tion) are similar if the original MB unfolding procedure 
described in Sec. 14. H is used. 

We assume a fixed value of axial mass for hydrogen 
Ma — 1030 MeV and we minimize distribution for an 
effective axial mass for carbon (M^^-^) using the data for 
the reconstructed energy distribution. As shown in 
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FIG. 3: [Color online] The ratio 77 (see Eq. I14p as a function 
of the reconstructed kinetic energy of all nucleons in a final 
state measured by MiniBooNE together with our best fits 
for two cases: the MEC contribution was not included (solid 
fine) and the MEC contribution was included (dashed line) 



this distribution is essentially insensitive to the strange 
quarks contribution, so here we assume g% — 0. Best 
fits for both cases are presented on Fig. [2] 

In the case without the MEC contribution the value 



M 



eff _ 



1312 ± 93 MeV 



(18) 



was obtained with Xmin/^OF = 26.1/50 (a confidence 
level (CL) 99.8%). This value is consistent with the MB 
result from Q (M^^ = 1390 ± 110 MeV). 
Inclusion of the MEC contribution leads to 



Ma = 1003 ± 87 MeV 



(19) 



with xl^in/DOF = 25.8/50 (CL 99.8%). This resuh is 
consistent with the world average axial mass value and 
confirms that the difference between recent and older 
axial mass measurements can be explained by taking 
into account the 2p — 2h constribution. 

To see how the choice of the electromagnetic form fac- 
tors parametrization affects the results of our analysis, 
we applied the form-factors corrected by the two-photon 
exchange [H,!!!]. We obtained = 1322 ±80 MeV 

with xlrrn/DOF = 25.8/50 (CL 99.8%) for the model 
without the MEC contribution and Ma = 1002±83 MeV 
with xl^^n|DOF = 25.4/50 (CL 99.9%) with the MEC 
included. This means that our results do not depend 
on the choice of parametrization of the electromagnetic 
form factors. 

Repeating the computations under the assumption 
that the value of the axial mass is the same for a 
free proton and bound nucleons, we obtained Ma = 
1032 ±93 MeV with xLn/^OF = 25.7/50 (CL 99.8%). 

Using the data for ratio (Eq. [Til) we examined the 
strange quark contribution to the NCEL cross section. 



We assumed a fixed value of axial mass: 1030 MeV for 
hydrogen, 1003 MeV for carbon when the MEC contri- 
bution is included and 1312 MeV if it is not. Best fits 
for both cases are shown in Fig. |3l 

We found the strange quark contribution to be 

5^ = -0.04 ±0.45 

with x'Lin/DOF = 27.8/29 (CL 52.9%), when the MEC 
contribution was not included. This results is consistent 
with values pubHshed by MiniBooNE Q and BNL E734 
0. 

MEC events contribute mostly to the denominator of 
the ratio rj. Also lower Ma, makes the ratio rj smaller. 
To compensate for those two effects a lower value of 5^ is 
required. Indeed, with the MEC contribution included, 
we got the value 



5a 



-0.92 ± 1.13 



with xlmin/DOF = 30.6/29 (CL 38.5%). The obtained 
value is still consistent with zero. 



6. CONCLUSIONS AND OUTLOOK 

The impact of the MEC contribution on the analysis 
of the MB data for neutrino NCEL scattering on CH2 
was investigated in detail. This is the first analysis of 
this kind yet. The final results confirm that recent large 
axial mass measurements can be explained by the two- 
body current contribution to the cross section. 

The value of Ma = 1003 ± 87 MeV was obtained for 
carbon nucleons after assuming Ma — 1030 MeV for free 
protons. The assumption that Ma is the same for car- 
bon and hydrogen led to value of Ma = 1032 ± 93 MeV. 

Our study relies extensively on the performance of 
the NuWro nucleon cascade model. The final results in- 
dicate that the model works well. This is encouraging, 
because an unambiguous separation of the two body cur- 
rent contribution from genuine CCQE events can only 
be done by looking at final state nucleons [H, |4^ . This 
is quite challenging, because a good control of nucleon 
final state interactions (FSI) effects is necessary. In 
particular, a good resolution of final nucleons with a 
low threshold for momenta of reconstructed tracks is re- 
quired. There is a hope that a liquid argon detector is 
able to provide necessary information and can be used 
to benchmark the performance of MC cascade models 
0. 

The strange quarks contribution to form factors was 
extracted using a sample of NCEL proton enriched 
events. The result is consistent with zero, however, due 
to limited accuracy of the data, the uncertainty of the 
best fit value of Qa is very large. 

Recently MiniBooNE made public preliminary results 
from the antineutrino NCEL analysis 18[. They are 
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supplementary to those discussed in this paper. When 
the data become available a combined analysis will have 
more potential to measure the MEC contribution and 
put more constraints on . 

We would like to notice that there is an interesting 
idea for an alternative measurement of the NCEL cross 
section described in 42|. The authors investigate the 



relation between the rate of observed 7 rays coming from 
nuclear deexcitation in water- Cherenkov detectors and 
the NCEL cross section. 
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